Superconductive properties of thin dirty SN bilayers 
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The theory of superconductivity in thin S'A'^ sandwiches (bilayers) in the diffusive limit is developed 
within the standard Usadel equation method, with particular emphasis on the case of very thin 
superconductive layers, ds <^ dN- The proximity effect in the system is governed by the interlayer 
interface resistance (per channel) pint- The case of relatively low resistance (which can still have 
large absolute values) can be completely studied analytically. The theory describing the bilayer in 
this limit is of BCS type but with the minigap (in the single-particle density of states) Eg <^ A 
substituting the order parameter A in the standard BCS relations; the original relations are thus 
severely violated. In the opposite limit of an opaque interface, the behavior of the system is in many 
respects close to the BCS predictions. Over the entire range of pint, the properties of the bilayer 
are found numerically. Finally, it is shown that the results obtained for the bilayer also apply to 
more complicated structures such as SNS and NSN trilayers, SNINS and NSISN systems, and SN 
superlattices. 

PACS numbers: 74.50.+r, 74.80.Fp, 74.20.Fg. 



I. INTRODUCTION 

It is well known that the majority of metallic super- 
conductors is well described by the classical BCS the- 
ory of supcrconductivity.EJ One of the main qualitative 
features of the BCS theory is a simple relation between 
the superconductive transition temperature Tc and the 
low-temperature value of the energy gap for s-wave su- 
perconductors: A(0) = 1.76 Tc. Experimentally, vio- 
lations of this simple relation are considered as a sign 
of some unusual pairing symmetry or even of a non- 
BCS pairing mechanism; many theories of unconven- 
tional superconductivity were developed during the last 
decade, mainly in relation with high-Tj, materials. Re- 
cently, an evident example of such a violation of the 
BCS theory nredictions was found in experiments by Ka- 
sumov et al.p who studied current- voltage characteristics 
of a carbon nanotube contact between two metallic bi- 
layers (sandwiches) made of ordinary metals, tantalum 
and gold. The observed value of the low-temperature 
Josephson critical current is 40 times larger thap the 
maximum expected (Ambegaokar-Baratoff) valueu I^. = 
7rA(0)/2ei?tubo, where the energy gap of the bilayer A(0) 
is estimated from its transition temperature. The source 
of such discrepancy is not clear at present. The most 
recent experimentsa demonstrate the existence of intrin- 
sic superconductivity in carbon nanotubes. However, 
the discrepancy could be due to unusual superconduc- 
tive properties of the bilayers. The aim of the present 
paper is to investigate these properties. 

An essential feature of the experiment [ ||] was that the 
superconductive layer in the bilayer was very thin {ds /d-N 
= 5 nm/100 nm = 1/20). In the present paper, we inves- 
tigate such a bilayer both analytically and numerically, 
calculating quantities characterizing the superconductiv- 
ity in this proximity system: the order parameter A, the 
density of the superconducting electrons n, the critical 



temperature Tc, the (mini)gap Eg in the single-particle 
density of states (DOS) and the DOS i^{E) itself, the 
critical magnetic field He parallel to the bilayer and the 
upper critical field Hc2 perpendicular to the bilayer. In 
our calculations, the parameter controlling the strength 
of the proximity effect is the (dimensionless) resistance 
of the SN interface per channel pint , which is related to 
the total interface resistance i?int as 



^int — 



RqPint 

2A^ch 



(1) 



where Rg ~ hje^ is the quantum resistance, and A^ch = 
^/(Af/2)^, with \p being the Fermi wave-length, is the 
number of channels in the interface of area A. We choose 
\f referring to the S layer. 

Our results show that the values of the interface re- 
sistance can be divided into three ranges: (a) at large 
resistance, many characteristics of the superconductor 
[A, n, Tc, v{E), He, Hc2] are almost unaffected by the 
presence of the normal layer, i.e., this is the BCS limit 
(however, we note that Eg does not coincide with the or- 
der parameter, and even vanishes as pint increases); (b) 
at low resistance, the theory describing the bilayer is of 
BCS type but with the order parameter A substituted 
by the minigap Eg (for instance. Eg = 1.76 Tc, whereas 
Eg <C A); the original BCS relations are thus severely 
violated; (c) at intermediate resistance, the behavior of 
the system interpolates between the above two regimes. 

The paper is organized as follows. In Section |l[ we for- 
mulate the standard technique of the Usadel equations 
for dirty systemsjj and introduce a convenient angular 
parameterization of the quasiclassical Green functions 
entering these equations. In Section HI, we apply the 



Usadel equations to the thin bilayer that we intend to 
discuss, and present numerical results for A, n, Tc, and 
i^{E). We start an analytical analysis of the Usadel equa- 
tions for the bilayer by calculating the minigap Eg in the 
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density of states in the two limiting cases of high and low 
interface resistance (Section Then, in Section 0, 

we elucidate the structure of the theory describing the 
system in the so-called Anderson limit (of relatively low 
resistance), finding A, n, Tc, Eg, and ^{E). The critical 
magnetic fields He (parallel to the bilayer) and Hc2 (per- 
pendicular to the bilayer) are ca lcula ted in Sections VI 
and VII, respectively. In Section VIII , we show that the 
results obtained for the bilayer also apply to more compli- 
cated structures such as SNS and NSN trilayers, SNINS 
and NSISN systems, and S'Af superlattices. The relation- 
ship between our results and the experiment [ |^ that 
has stimulated our research is discussed in Section IX. 
Finally, we present our conclusions (Section ^). 



II. METHOD 

A. Usadel equation 

Equilibrium properties of dirty systems are described^ 
by the quasiclassical retarded Green function R{r,E), 
which is a 2 X 2 matrix in the Nambu space satisfying 
the normalization condition — 1. The retarded Green 
function obeys the Usadel equation 



DV{RVR) + i[H, R] = 0. 



(2) 



Here the square brackets denote the commutator, D ~ 
vl/3 is the diffusion constant with v and I being the Fermi 
velocity and the elastic mean free path, H = Ea^ + A(r) 
with E being the energy, whereas cTz (the Pauli matrix) 
and A(r) are given by 



1 

-1 



A = 



A 

-A* 



(3) 



The order parameter A(r) must be determined self- 
consistently from the equation 



A(r) 



dE tanh. 

2T 



R{v,E) - R{v,-E) 



o.d. 

(4) 



where the subscript o.d. denotes the off-diagonal part, 
vq = m^v/TT^ is the normal- metal density of states at the 
Fermi level, A is the effective constant of electron-electron 
interaction in the S layer (whereas we assume A = and 
hence A = in the Allayer), and integration is cut off at 
the Debye energy lo£) of the S material. 

Equation (^ should be supplemented with the appro- 
priate boundary conditions at an interface, which readQ 



(^RlVnRl) = <Jr (Rr^nRr) = ^ 



Rl, Rr 



(5) 



where the subscripts I and r designate the left and right 
electrode, respectively; a is the conductivity of a metal in 



the normal state, and g-„it = Gint/A (with Gint = l/i?int) 
is the conductance of the interface per unit area when 
both left and right electrodes are in the normal state. 
Vn denotes the projection of the gradient upon the unit 
vector n normal to the interface. 

The system of units in which the Planck constant and 
the speed of light equal unity {ti = c = 1) is used through- 
out the paper. 



B. Angular parameterization of the Green function 

The normalization condition allows the angular param- 
eterization of the retarded Green function: 



ie'f sin ( 



— cost 



(6) 



where 6 = 6{y, E) is a complex angle which characterizes 
the pairing, and Lp — (/j(r, E) is the real superconduct- 
ing phase. The off-diagonal elements of the matrix R 
describes the superconductive correlations, vanishing in 
the bulk of a normal metal {9 — Q). 
The Usadel equation takes the form 



D 



D , ,2 

iE — — (V(^) cos( 



sin 6* -I- lAI cos 6* = 0, 



(7a) 



V (sin2 6'V(^) = 0. (7b) 
The corresponding boundary conditions are 

O'l'^n^l = dint [C0s{(pr — ifil) COS 01 SlllOr — SmOl COS 0^] , 

(8a) 

CrVn^r = ffint [cOsOl SluOr — COs{ipr — fl) siu 0; COS 0^] , 

(8b) 



Ul SIT? OlVafl — (Tr siu^ OrVnfr 

= gint sm{ipr - (fi) sin 01 sinBr- (8c) 

The self-consistency equation for the order parameter 
A(r) takes the form 



A = 



I'oA /■'^^ f E\ 



(9) 



The above equations are written in the absence of an 
external magnetic field. To take account of the magnetic 
field, it is sufficient to substitute the superconducting 
phase gradient in the Usadel equations (0) by its gauge 
invariant form 2mv = Vyj -I- 2eA, where A is the vector 
potential and v denotes the supercurrent velocity. 

Physical properties of the system can be expressed in 
terms of the pairing angle 9{r,E). The single-particle 
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density of states i/(r, E) and the density of the supercon- 
ducting electrons n(r) are given by 



III. USADEL EQUATIONS FOR A THIN 
BILAYER 



= f Re [cos 6] , 
2mtT 



J dE tanh Im [sin^ 9] , 



(10) 
(11) 



where m and e are the electron's mass and the absolute 
value of its charge. The total number of single-particle 
states in a metal is the same in the superconducting and 
normal states, which is expressed by the constraint 



dE[i^{r,E)-i^o]=0- (12) 



C. Simple example: the BCS case 

The simplest illustration for the above technique is the 
BCS case, when the order parameter A(r) = ^bcs is 
spatially constant. Its phase can be set equal to zero, 
ip = 0. Then the Usadel equations (0) are trivially solved, 
and we can write the answer in terms of the sine and the 
cosine of the pairing angle: 



sine BCs{E) = 



cos 0BCsiE) = 



iAscs 



^E^ - A2 



BCS 



E 



(13a) 
(13b) 



An infinitesimal term iO should be added to the energy 
E to take the retarded nature of the Green function R 
into account, which yields 

Im [sin^ OBcsiE)] = ^Abcs S{E - Abcs)- (14) 

The usual BCS relations are straightforwardly ob- 
tained from Eqs. (^, (|l^), ( pT| ) (for simplicity, we con- 
sider the case of zero temperature): 



(15) 

if £; < Abcs 

if E> Abcs ' ^^"^^ 



Abcs = 2cj_d exp ( ) , 



VBCs{E) 



0, 
ma 



E 



BCS 



nBCS = tt^-Abcs- 



(17) 



The critical temperature must be determined from 
Eq. (0) with vanishing A(Tc); the result is 



BCS 



Abcs 

where 7 ~ 1.78 is Eulcr's constant 



(18) 



Let us consider a SN bilayer consisting of a normal 
metal {—dN < z < 0) in contact (at z = 0) with a super- 
conductor (0 < z < ds). We assume that the layers are 
thin (this assumption will be discussed in Sec. IX) and 
can be regarded as uniform, which allows us to set the 
order parameter A equal to a constant in the supercon- 
ductive layer (we choose its phase (p equal to zero). At 
the same time, we suppose that electron-electron inter- 
action is absent in the normal layer: A = 0, hence A = 0, 
although the superconductive correlations {9 ^ 0) exist 
in the layer due to the proximity effect. The Usadel 
equations (0) take the form 



Dn d^9N 
2 dz^ 



iEsm9N — 0, 



Ds d^Os 
2 dz^ 



iE sin 9s + A cos 9s = 0, 



(19a) 
(19b) 



where 9n and 9s denote the pairing angle 9 at z < and 
z > 0, respectively. 

The boundary conditions (0) reduce to 



d9N d9s . , , 

CTjv-^ = o-S^— = ffint sm(6'5 - 9n). 
az az 

Equations (|l9|) can be integrated once, yielding 



(20) 



Dn (d9 



N 



4, \ dz 



iEcos9M — /jv, (21a) 



D / f)9 \ 

~f[~d7j -*^^os0s + Asin05 = (21b) 

The functions fN{E) and fs{E) are determined from the 
boundary condition 39 /dz = at the nontransparent 
outer surfaces of the bilayer, which give 

fNiE) = ~iEcos9Ni-dN,E), 

fs{E) ^ ^iEcos9s{ds,E) + Asin9s{ds,E). (22) 

Let us denote 9n{E) = 6'jv(-dw, -B), 9s{E) = 9s{ds, E). 
Because of the uniformity of the layers, the functions 
6n{z,E) and 9s{z,E) are nearly spatially constant. 
However, in order to determine them, we should take 
account of their weak spatial dependence and make use 
of the boundary conditions at the SN interface. Substi- 
tuting 



9N{z,E)^9NiE)+S9N{z,E), 
9s{z,E)^9s{E)+S9s{z,E) 



(23) 

into Eqs. ( pT| ) and linearizing them with respect to 
\69Niz,E)\, \69siz,E)\ < 1, we find the solution. Fi- 
nally, boundary conditions at the SN interface lead to 

- iTNEsm9N{E) = iTsEsin9s{E) + tsAcos9s{E) 
^sm[9siE)-9NiE)], (24) 
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where we have denoted tn = '^aNd-N / gint, ts = 
'^.CTsds I Ds 5int • Using the definition of the interface resis- 
tance per channel (|l|), we can represent these quantities 
as 



TN 



^TT 2 Ant' 



TS = 27r— Pint, (25) 



with vn and vg being the Fermi velocities in the N and 
S layers. The ratio /ts = ynAn l^sds, which is inde- 
pendent of interface properties, can also be interpreted 
as the ratio of the global densities of states (per energy 
interval) in the two layers considered. 



TN 
TS 



AdNVON 

Adsvos ' 



(26) 



The latter interpretation will prove useful for further 
analysis. 

Having solved the boundary conditions (p^ ) , we can de- 
termine all equilibrium properties of the system (because 
knowledge of 9s, On implies knowledge of the retarded 
Green function R). 

A useful representation of the boundary condi- 
tions (041) is obtained as follows. Excluding 9n{E) from 
Eq. (|2J)7 we arrive at a single equation for the function 
9s{E), which can be written, in terms of Z = exp{i9s), 
as a polynomial equation 

iCeZ^ + C'sZ^ + iCiZ"^ + C^Z^ + iCa^^ + CiZ + iCo 
= (27) 

with real coefficients 

\TN 

C5=[l- [tnE]- 
Ca = -TnE 



2 r 



1 + 



TS 
TN 



1 + — 



TS 
TN 



C. = 2 - 2 



1 - [tnEY 



TS 
TN 



A 



C2 

Co 



-TnE 



Ts 
TN 



-TnE 



tnE) 

Ts 



TN 



I TS 

\tn 

2 r 

1 - 



A 

2 

A' 



e 



A 

E 



(28) 



During further analysis, the choice between the bound- 
ary conditions in the forms ( p4[ ) and (^^ will be a matter 
of convenience. 



A. Critical temperature 

The critical temperature of the bilayer Tc is defined 
from the condition of vanishing of the order parameter 
A. Near Tc, the superconducting correlations are very 
small, \9\ <C 1; nevertheless, the self-consistency equa- 
tion (||) has a nonzero solution A ^ 0. 

Linearizing the boundary conditions ( p4[ ) with respect 
to 9n and 9s, we readily find the solution: 



ME) = 4(1- 

E \ Ts 



tn 



TN - itstnE 



(29) 



Substituting this into the self-consistency equation (||) 
and simplifying both its sides by A, we obtain an equa- 
tion determining T^, which can be cast into the form 



In- 



j^BCS 



TN 



Tc 



TS + TN 



"^'2 2ttTcTstnJ 




, (30) 



where ipix) denotes the digamma function. A similar 
formula (except the logpjithmic term in the r.h.s.) was 
obtained by McMillanB^lia (see also Ref. [|l|]). The loga- 
rithmic term in the r.h.s. takes account of the finiteness 
of the Debye energy wd; it becomes important only in the 
limit of a perfect interface (the Cooper limit), i.e., when 
TsTNi^o/iTs + Tn) <C 1. Equation ( |30| ) can be solve d 
numerically over the entire range of pint (see Sec. Ill B| ) ; 
the analytical solution can be found in limiting cases (see 
Sec.0). 



B. Numerical results 

The solution of Eq. (^^ can be found numerically. To 
this end, we solve the system of two nonlinear equations 
for the functions IleZ{E) and lmZ{E), using the modi- 
fied Newton method with normalization. 

The solution depends on the bilayer's parameters: the 
thicknesses of the layers, characteristics of materials con- 
stituting the bilayer, and the quality of the 5A'^ interface. 
This dependence enters Eqs. (|27|), (^ ) via tn and ts- 
For numerical calculations, we assume the characteris- 
tics of the bilayer ip be the same as in the experiment 
by Kasumov et aln The superconductive layer is made 
of tantalum, ds = 5 nm, and the normal layer is made 
of gold, dN = 100 nm. Approximate experimental val- 
ues of the conductivities arellj as = 0.01 /ifi"^ cm~^ and 
(Tat = 1 /ifJ^^ cm^^. In order to calculate the Fermi char- 
acteristics of tantalum and gold, we use the values of the 
Fermi energy EpiTa) =r-k^ eV, i;^F(Au) = 5.5 eV, and 
the free electrons model.Ej 
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FIG. 1. Dependence of the order parameter in the S layer 
A, of the superconducting electrons' density in the S layer 
ns, and of the bilayer's critical temperature Tc on the inter- 
face resistance per channel pint at zero temperature. All the 
quantities are normalized by the corresponding BCS values. 
The discrepancy between the curves implies a violation of the 
BCS relations between A, ns, and Tc- 
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FIG. 2. Zoomed part of Fig. |l[ In the shown range of rel- 
atively small resistance pint, the BCS relations between A, 
ns, and Tc are severely violated. The upper and lower graphs 
differ only in the scaling of the ordinate axis (normal and 
logarithmic, respectively). 
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FIG. 3. Temperature dependence of A and ns at Pint — 
90 and 110. The temperature is normalized by the critical 
value Tc, which depends on pint; A and ns are normalized 
by their zero-temperature values. For comparison, the same 
dependence is also plotted for the BCS case. 



:80 







o 

> 


2.0 


0.5 


A/A^cs 








v(E)/ 


1.5 


0.4 


ng/ns'^'' 










1.0 


0.3 


J / J BCS / 
c ' c 










0.5 


0.2 












0.0 


0.1 ■ 












3.0 


0.0 - 












2.5 


8 


90 100 110 120 


130 


140 


150 


o 

> 


2.0 



S and N layers 



: 110 



S layer 

N layer 



0.5 1.0 1.5 2.0 



0.0 0.5 1.0 1.5 



Pin. = 150 




- 1 
1 

1 1 


^^m = 2000 _ 


1 7 
1 

- 1 / 
1 J 
_ 1 


S layer 

N layer 



0.5 1.0 



E/A E/A 
FIG. 4. Energy dependence of the single-particle DOS, nor- 
malized by the normal-metal DOS, v(E)Ivq, in the S and 
layers at pint = 80, 110, 150, and 2000. The energy E is nor- 
malized by the order parameter A, which is different in all the 
four cases. [To avoid confusion, we note that actual (in ab- 
solute units) relation between the minigaps in the DOS does 
not correspond to what is seen from the Figure.] The Figure 
demonstrates drastic difference of the DOS in the bilayer from 
the BCS case [Eq. (Fol)]. 
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Once the parameters have been specified, the solution 
of Eq. (^ ) depends only on the interface resistance pint- 
Having found the function Z{E) [which is equivalent to 
finding 6s{E)], we start from the case of zero temper- 
ature, T = 0, and study the dependence of the order 
parameter A and of the superconducting electrons' den- 
sity in the 5* layer ^5 [Eqs. (||), (|ll|)] on pint- The results 
are plotted in Fig. |l|, where we also show the dependence 
of the critical temperature Tc, determined from Eq. (30), 

on Pint- 

The suppression of A, 715, and T^, in comparison to 
their BCS values in the 5 layer, is a natural consequence 
of proximity to the normal metal. At the same time, 
there is a possibility of BCS-like behavior, which implies 
the BCS relations between the suppressed quantities and 
the coincidence of the three curves plotted in Fig. ^ How- 
ever, the curves split, and the difference between them 
is largest for relatively small values of pint- Figure ^ 
presents the range 80 < pint < 150 on a larger scale. 

Figure || shows the temperature dependence of the or- 
der parameter A and of the superconducting electrons' 
density in the S layer ^5. Although the smaller pi„t 
the further it is from the BCS limit (corresponding to 
Pint 00), we observe that at pint = 90 the curves are 
closer to the BCS behavior than at pint = HO. An ex- 
planation of this feature is given in Sec. 0. 

Finally, the energy dependence of the single-particle 
density of states in the S and layers vs,n{E) is plotted 
in Fig. 1^. The density of states in the bilayer is qualita- 
tively different from the BCS result (p^. In particular, 
there is a minigap Eg in the density of states at energies 
much smaller then Aspg, and even much smaller then A 
in the bilayer. In the next Section, we find this minigap 
analytically in the two limiting cases of small and large 
interface resistance pint- Another feature which can be 
seen from Fig. ^ is that the order parameter A plays the 
role of a characteristic energy scale of the system only in 
the limit of large pint (see the case pi„t = 2000) . Other- 
wise, no peculiarity in the DOS is observed at i? = A. 
Some other aspects of the DOS behavior will be discussed 



in Sees. IV 







IV. MINIGAP IN THE DENSITY OF STATES 

In principle, the energy dependence of the single- 
particle density of states is different in the Sand A^layers. 
At the same time, the gap in the DOS is a property of 
the bilayer as a whole; the gap is spatially independent 
because there is no localization in the system and each 
electronic state extends over the whole bilayer. 

The presence of the gap thus implies that the density of 
states vanishes in both layers, when the energy is below 
the gap: 

Re [cos 6*5] = Re [cos6'jv] = 0, (31) 

leading to 6*5 = ^ -f ids, Sn — ^ + ii}N, with real -ds and 
■&N- In this case, Eqs. ( p^ ) can be written as 



tanhi? 



N 



sinhi^S -I- TpfE 
cosh ??5 ' 



cosh^Jjv — — — coshiJs + — ^ 

TN TnE 



sinh . 



(32a) 
(32b) 



Assuming that sinh 3> TjyE at small energies, from 
Eqs. (§^) we obtain = ^s, and, finally. 



cos 6*5 — cos Sat = 



E 



E'^^El 



(33) 



with Eg = tsA/(t5 -I- tn). This is a BCS-hke result 
[cf. Eq. ( |l3b| )], although the order parameter Abcs is 
substituted by the minigap Eg. The assumption is read- 
ily checked, and we obtain 



Eg = ^^A, if ™^ « 1. (34) 

Tn +TS Ts + Tn 



Now we proceed to the opposite limit of large interface 
resistance. Assuming sinh-iJs <C tnE and sinh ^ 1, 
we solve Eqs. ( ^ ) and finally obtain 



cos Us = 



E .E 
rsA^E^ - \It% " * A' 
E 



cos 



- 1/r 



(35a) 
(35b) 



The assumption is readily checked, and the result is 



Eg = — , if 

TN 



tstnA 

Ts + TN 



> 1. 



(36) 



Equations ( |34D and (|3g) imply that Eg is a nonmono- 
tonic function of the interface resistance: with increase 
of Pint , it first increases at small pint [Eq. (|3j)] and then 
decreases at large pint [Eq. (p6|)]. Therefore, Eg reaches 
its maximum at some intermediate value of pint, corre- 
sponding to tstn A./ {ts + Tjv) ~ 1, hence pint ~ 140. 
Numerical results for Eg are shown in Fig. ^. 

At first sight, vanishing of the minigap Eg in the limit 
of an opaque interface seems to contradict the general 
tendency to the BCS behavior. However, this contradic- 
tion is more apparent than real. Actually, the DOS curve 
for the S layer does approach the BCS result (|l6|) in this 
limit, showing the standard peculiarity at ~ Abcs- 
At the same time, below Abcs, the DOS curve sharply 
drops to very small values (which are still finite in con- 
trast to the BCS case), and turns to zero only at E = Eg. 
Simultaneously, the DOS in the N layer approaches the 
(constant) normal-metal value. The tendency to such 
behavior is illustrated by Fig. B, the case pint — 2000. 
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100 150 200 250 300 350 400 

Pint 

FIG. 5. Minigap in the single-particle density of states Eg 
versus pint. The minigap is normalized by the BCS gap value 
Abcs- Eg is a nonmonotonic function of pint, reaching its 
maximum at pint ~ 160. The inset shows -Eg (pint) on a wider 
(logarithmic) scale over Pint . 

The results (^^, ( ^ ) are valid not only below the mini- 
gap but also right above it (in which case the real parts 
of 9s and On differ from 7r/2), providing a comparison 
between the DOS in the S and TV layers. Equation (|3^ ) 
demonstrates the equality of the DOS in the two layers 
at relatively small pint (see Fig. U, the case pint = 80). 
Proceeding to the limit of an opaque interface, we should 
note that the approximation which lead to Eqs. ( |35| ) fails 
in a narrow vicinity of Eg [this fact does not affect the 
result for the minigap ( pq ) itself, but the incorrect diver- 
gence of the DOS ai E = Eg disappears]. Equations ( |35| ) 
demonstrate that outside this region, at energies right 
above the minigap, the DOS in the layer exceeds the 
DOS in the S layer by the large factor r^A (see Fig. ^, 
the case pint — 2000). 

Finally, we note that the results of the preappt Section 
are similar to those obtained by McMillanoEl (see also 
Ref. [Inl). 



V. ANDERSON LIMIT 

In the limit of relatively low interface resistance (the 
Anderson limit), the theory describing the bilayer can be 
developed analytically. The condition defining this limit 
is rsA, tatA < 1. 

First of all, we need to determine 9{E) [or Z{E)] solv- 
ing Eq. (|2^) [or Eq. (^7|)] over the entire range of energies 
E. 

In the region > A, the solution of Eq. (27) can be 
written as Z — 1 + dZ, with \SZ\ <^ 1. Keeping terms up 
to the second order in SZ , we obtain 



6Z^ 



A(l - iTNE) 



itnE^ 



(37) 



This result is general in the sense that it is valid for ar- 
bitrary values of pint . 

At i? < A, the same calculation as for the minigap 
leads to the result 



iEa 



E^-E^ 



(38) 



with the minigap Eg given by Eq. (|34|). [To avoid con- 
fusion, we note that under the less strict limitations for 
TsA, tnA used in Eq. (H), the BCS-like results (|3|) 
and (^ are valid only up to energies of the order of Eg.] 
Now Im[sin05] is readily calculated, and in the case 
of zero temperature, T = 0, the self-consistency equa- 
tion (||) can be solved, yielding 



A 



A 



BCS 



TS 



2w_D y \Ts + Tn ) 



tn /ts 



(39) 



The relation between the order parameter A and the 
minigap Eg is given by Eq. (|3^), which immediately 
yields 



E„ 



A 



BCS 




tn/ts 



(40) 



In the limit of a perfect interface (the Cooper limit), 
which is defined by the condition tstn^jJo I (ts^t^n) <C 1, 
Eq. (^) reproduces classical Cooper's resultliJ general- 
ized to the caSjC of different Fermi parameters in the S 
and layers:! 



Egipint 0) = 2tjDexp 



with the effective interaction parameter 



TS 



TS + tn 



(41) 



(42) 



This parameter can be considered as a result of averaging 
with the weighting factors r, which are proportional to 
the global normal-metal DOS per unit interval of energy, 
T oc Advo (note that the interaction parameter is zero in 
the Allayer). 

At the same time, we would like to emphasize that 
the Anderson limit docs not reduce to the Cooper limit 
with small corrections. On the contrary, due to the re- 
lation A <C cj_D, the Cooper limit's condition is not sat- 
isfied over the most part of the Anderson limit's valid- 
ity range; therefore, the minigap Eg and the quantities 
calculated below differ drastically from the Cooper limit 
expressions. 

Now we proceed to calculate the density of the super- 
conducting electrons in the S layer ng. On this way, we 
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immediately encounter the problem that the above solu- 
tion (^8|) of Eq. at i? < A is not accurate enough for 
our purpose. In fact, as we will see below, the principal 
contribution to the integral ([III ) determining ns comes 
from a narrow region of energies near Eg. At the same 
time, Eq. (|38| ) yields Im [sin^ 9s] — and hence no con- 
tribution at all from E < A. We thus have to calculate 
a correction to Eq. (js^). Assuming this correction to be 
small, we linearize Eq. and obtain 



iEr, 



with 



sm 



X{E) = 



E^^E^ 



X{E), 



tnE„ 



(43) 



(44) 



which is valid for all < A except for a narrow vicinity 

The accurate consideration of the minigap's vicinity is 
possible due to the fact that \Z\ <C 1 in this region. We 
define a new dimensionless quantity e as 



E 



(rsA) 



2/3 



2 1 



tn 



4/3 



(45) 



and consider the region \E — Eg\ <C Eg. The function 
Im[sin'^05], which determines ns, has a peak at e ~ 1. 
An analysis of the coefficients (^8|) shows that only the 
terms that contain zeroth ^ fi rst, and third order in Z 
should be retained in Eq. ( |27| ). Then, after rescaling 



Z 



(rsA) 



1/3 



1 



2/3 



(46) 



we obtain a cubic equation for the function Y{e): 

-eY + 0, (47) 

which can be solved analytically. 

The density of the superconducting electrons at 
zero temperature is now readily calculated: 



ns 

^BCS 



E„ 



11 (rsA)' 



5/6 



^BCS 



2-K 



4/3 



Its A 



TT ( 1 + ^ 



1 J- "^s 

I^ + 21ni±^ 

Ts TsA 



(48) 



with Eg given by Eq. (|40|). The first term in the curly 
brackets is the principal one; the two other terms be- 
come comparable to unity only near the upper limit of 
applicability of Eq. (El). 



The critical temperature of the bilayer Tc is determined 
byEq. (^o|). In the Anderson limit, {ts+tn)/tstn ^ Tc, 
and, with the use of the asymptotic form of the digamma 
function ip{^) ^ Inx at x 3> 1, we obtain 



E„ 



TBCS A 



BCS 



(49) 



with Eg given by Eq. (j^j. Interestingly, this result ex- 
plains a discrepancy in the formulas-, for Tc of a thiu 
bilayer that were found by CooperliJ and McMillan. □ 
This discrjepancy is discussed in the classical paper by 
McMillan JS'lI We conclude that both cited results are cor- 
rect, but their applicability ranges are different, although 
within the Anderson limit. Cooper's result corresponds 
to a perfect interface, tstn^d / {ts + ta?) <C 1, whereas 
McMillan's formula applies in the case tstmojd I {ts + 
tat) » 1. 

Now we can discuss the general structure of the the- 
ory describing the bilayer in the Anderson limit. In the 
limit Pint — ^ 0, our results for the pairing angle 9 (which 
is constant over the entire bilayer, 9 = 9s = 9^) yield 
expressions which can be obtained from the BCS ones 
[Eqs. dl^), (p^J if we substitute the BCS order param- 
eter Abcs by the bilayer's minigap Eg. At pint > 0, 
corrections to this simple result are small while the An- 
derson limit's conditions are satisfied. Therefore, we ob- 
tain a BCS-type theory with Eg substituting Abcs hi 
all formulas. 

The results of this Section immediately explain the nu- 
merical results in the limit of relatively small pint , shown 
in Fig. 1^. As we have found, the Anderson limit implies 
the following relations between the quantities under dis- 



Er, 



-Tc 



mas 
ns = T^—^Eg, 



(50) 
(51) 



which substitute Eqs. (|l|) and (|T^). For the Ta/Au bi- 
layer to which the numerical results refer, the Anderson 
limit is valid at pint < 80 (we see that the values of pint 
can be large although they are relatively small). There- 
fore, approaching pint = 80, the curves ns/ng^^ and 
Tc/T^'-^^ tend to coincide, and A/ Abcs exceeds them 
by the large factor (1 + /ts) ~ 15. 

The temperature dependence of A and ns, shown in 
Fig. ^, is quite different at pint = 90 and pint = 110; at 
Pint = 90, the curves are much closer to the BCS behav- 
ior. This is also explained by approaching the Anderson 
limit, where the curves coincide with the BCS ones. 

The DOS in the S and N layers coincide in the Ander- 
son limit [Eqs. (^Sf), (^8|)], showing the standard BCS-like 
peculiarity at E — Eg. The tendency to such behavior is 
illustrated by the DOS curve for pint = 80 in Fig. 0. 
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VI. PARALLEL CRITICAL FIELD 

We proceed to calculate the critical magnetic field H, 
directed along the plane of the bilayer. As it was men 
tioned in Sec 



II B 



in the presence of an external mag- 
netic field, the superconducting phase gradient in the 
Usadel equations (0) must be substituted by its gauge 
invariant form, which can be expressed via the super- 
current velocity v. The spatial distribution of v in the 
bilayer can be found as follows. 

Let us direct the x-axis along the magnetic field H. 
The supercurrents j = — env are directed along the bi- 
layer and perpendicularly to H, i.e., j = (0,j(z),0) and 
V = (0, 0). Near He, the magnetic field inside the bi- 
layer is uniform, so the vector potential can be chosen as 
A = (0, — zi?, 0). The supercurrent velocity distribution 
is determined by the equation V x v = eH/m. Another 
essential point is the continuity of v at the SN interface, 
which follows from the continuity of the superconducting 
phase ip [see the boundary condition (pel)]. The result is 



eH 

v[z) =Vo z, 

m 



(52) 



where is the supercurrent velocity at the interface, 
which must be determined from the condition that the 
total charge transfer across the bilayer's cross-section is 
zero: 



j{z)dz = 0, 



— dp. 



leading to 



N 



2m / nsds + n^d 



N 



(53) 



(54) 



The density of the superconducting electrons is constant 
in each layer [ns and rtAr). 

Near He, the superconducting correlations are small, 
\6\ <C 1, and the Usadel equation (|7a| ) for the paring 
angle 6{z, E) can be linearized: 



N 



2 dz^ 



UE -2m^DN^'^) 9n = 0, 



Ds d^Os 
2 dz^ 



{iE - 2m^Ds v^) Bs + | A| = 0. 



(55a) 
(55b) 



At the same time, the second Usadel equation (7b) is 
trivial: its l.h.s. is proportional to 



V (sin^ Ov)^ sin 29 V9 v + sin^ 9 Vv, 



(56) 



where both terms vanish due to the fact that V0 is di- 
rected along the z-axis whereas v is parallel to the y-axis. 

The pairing angle 9 is almost spatially constant in each 
layer; this allows us to average each of Eqs. (55) over the 
thickness of the corresponding layer, obtaining 



o9n 




2dN 


dz 


z=0 


Dn 


d9^ 




2ds 


dz 


z=0 


Ds 



{EN~iE) 9n, 



|A|], 



where 



En = 2m^DN (v2(z)) 
Es^2m^Ds{w^{z)), 



(57a) 
(57b) 

(58) 



are i?-dependent energies. Using Eqs. (|52|), (|54|), we ex- 
press them via He and the densities of the superconduct- 
ing electrons: 



Es 



6 



dl + 3d^ 



n% jds + d^) 
{nsds + UNd^f 



, (59) 



and En is obtained by the interchange of all the 5* and 
indices. 

Substituting ( [S^ ) into the boundary conditions ( [20| ) 
(which should be linearized), we find 

9n = Ts |A| / {tsEs + tnEn + tstnEsEn - tstnE"^ 

-iE[Ts + TN + TsTN{Es + EN)]}, (60a) 
9s = {I + tnEn -iTNE)9N. (60b) 

The order parameter A cancels out from the self- 
consistency equation (^. However, the resulting equa- 
tion alone does not suffice for determining Hc{T) because 
it contains Es and En, which are functions of un /ns. 
Therefore, to obtain a closed system, we must consider 
the self-consistency equation together with the equation 
determining the ratio njv/«s; the latter equation is ob- 
tained from Eq. ([Tl]). The resulting system of two non- 
linear equations for the quantities He and nN /ns is 



In- 



2uj 



D 



dE tanh — -r-TT^. 
|A| ■ 



^N /o°° dE tanh(^) Im^^ 



ns 



crsJ^dE tanh(^) lm9l 



(61a) 
(61b) 



with 9n and 9s given by Eqs. (|60|). The first equation of 
the system, Eq. (ala), can be written via the digamma 
functions, thus taking exactly the same form as Eq. (|79| ) 
below (which determines the perpendicular upper critical 
field) if we denote £s ^ Es + 1/ts, £n = En + 1/tn. 

In the limit pint oo, Eqs. ( |6l| ) lead to the BCS re- 
sult. In this case, the layers uncouple, the density of 
the superconducting elect rons in the N layer vanishes, 
nN/ns — *■ 0, and Eq. (61a) finally yields 



In- 



rpBCS 



T 




Ds [eHf^^'ds]' 
12ttT 



(62) 



which determines the parallel critical field He'-^^{T) of 
a thin superconducting film. Another immediate conse- 
quence of Eqs. (|6|) is the critical temperature of the bi- 
layer Tc, which can be found from the condition Hc(Tc) = 
0: in this case, Eqs. (|6l|) reproduce Eq. (30). 
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FIG. 6. Parallel critical field He, normalized by the BCS 
value, versus pint at zero temperature. The upper and lower 
graphs differ only in the scaling of the ordinate axis (normal 
and logarithmic, respectively). The nature of the steep be- 
havior of He at Pint = 120-123, which is best seen from the 
lower graph, is explained in the text. The inset shows Hc(pi-nt) 
on a wider scale over pint. 



o 

0.6 

X 




P:„t=121.2 

BCS curve 



T/T, 

FIG. 7. Temperature dependence of the parallel critical 
field He at pint = 121.2. The experimental value of Hc{0), 
analyzed with the use of the results shown in Fig. ^, sug- 
gests that this value of pint corresponds to the experiment 
by Kasumov et ala The critical field is normalized by its 
zero-temperature value, and the temperature is normalized 
by the corresponding Tc- For comparison, the same depen- 
dence is plotted for the BCS case. 



FIG. 8. Position of the stationary point zo of the super- 
current distribution versus pint at zero temperature. The co- 
ordinate Zo is normalized by the S layer thickness ds- The 
fast shift in zq from the center of the S layer at large pint to 
(nearly) the center of the A'' layer at small pint corresponds to 
the steep drop in He, shown in Fig. |^. 

The system of equations ( |6T| ) can be solved numeri- 
cally at arbitrary values of the temperature T and the 
interface resistance p-mt', the results for He are presented 
in Figs. I, 0. 

A remarkable feature of the function -ffc(pint) at zero 
temperature (Fig. ^) is the steep behavior of He at 
Pint = 120-123. This feature is due to rearrangement 
of the supercurrents inside the bilayer, which occurs in 
the following way. The supercurrent velocity changes 
across the thickness of the bilayer according to the simple 
linear law (|5^). This supercurrent distribution may be 
characterized by the position of the stationary point zq, 
where the supercurrent velocity is zero: v{zo) = 0, hence 
zq = mva/eH. At large values of the interface resistance 
Pint, the density of the superconducting electrons in the 
TV layer is very small, hn/us ^ 1, and the supercur- 
rents circulate only in the S part of the system; this case 
corresponds to 



Zo 



2 ■ 



(63) 



Then, while decreasing pint, a shift in zq occurs. Now 
the supercurrents in the S layer are not compensated (in 
the sense of the charge transfer); therefore, they must be 
compensated by the supercurrents in the A^ layer, which 
are enhanced due to significant increase in n^. This 
situation corresponds to the beginning of the drop in 
He- The ratio of the superconducting electrons' densi- 
ties grows rapidly, approaching the Anderson limit value 
n-N/n-s — o'n/o's (see Sec. VIA below); simultaneously, 
zq tends to 



asdl 



N 



2 [asds + OMdN) ' 



(64) 



10 



and the steep drop in He finishes. For the bilayer to which 
the numerical results refer, ds <C cIn and as <C ctn, so 
Eq. ( |6^ ) yields zq ~ — c?Ar/2. 

This scenario is illustrated by Fig. ^, which has been 
obtained numerically. 

The analytical solution of Eqs. (^) at zero tempera- 
ture in the Anderson limit is presented below. 



A. He at zero temperature in the Anderson limit 

In the zero-temperature Anderson limit (defined by the 
conditions tsEs, tj^En <SC 1) , th e ratio of the supercon- 
ducting electrons' densities ( |61b| ) becomes independent 
of the magnetic field, hn /us = ctn/cfs, and the self- 
consistency equation (61a) yields 



TsEs + tnEn 

TS + TN 



^BCS 



A 



BCS 



2uJD 



(65) 

which determines He- This result can be compared to 
the BCS case, which corresponds to the limit pint oo. 
In this case, the density of the superconducting elec- 
trons in the A^ layer vanishes, n^/ns 0, and the self- 
consistency equation yields 



E 



BCS 



A 



BCS 



(66) 



where Eg^^ is given by Eq. (|59| ) with tiat = 0. Finally, 



H. 



BCS 



S,BCS 



Ds 



^BCS 



(67) 



where $o = ^r/e is the flux quantum, and ^bcs is the 
correlation length in the dirty limit. 

Remarking that the r.h.s. of Eq. (pa l i s identical to 
Eg/ 2 with the minigap Eg given by Eq. (^(]|), we see that 
equation (|65|), determining the parallel critical field of 
the bilayer in the Anderson limit, is obtained from the 
BCS equation ( |6^ ) if we substitute the order parameter 
^Bcs by the minigap Eg (in accordance with the results 
of Sec. 0) and the i/-dependent energy Es^^ by the cor- 
responding averaged quantity (tsEs + t^En)/ {ts + tn)- 

The explicit result for the parallel critical field of the 
bilayer, obtained from Eq. (65), can be cast into a BCS- 
like form: 



He 



(68) 



The bilayer's correlation length ^ is the characteristic 
space scale on which the order parameter (or the pairing 
angle 9, or the Green function) varies in the absence of 
the magnetic field. In the Anderson limit (under discus- 
sion), the explicit formula for ^ is a natural generalization 



of the BCS expression [see Eq. (^)] which implies that 
Ds must be substituted by the averaged diffusion con- 
stant (D) and Abcs must be substituted (in accordance 
with the results of Sec. 0) by the bilayer's characteristic 
energy scale, the minigap Eg [Eq. (^)]: 




TsDs + tnD 



N 



Eg ' " ' TS +TN 



(69) 



The effective thickness of the bilayer in Eq. (^) is 



dcs — 



{<7sds + <JNdN) (o-s^l + o-jvd; 



N 



-SasCTNdsdN [ds + d 



N) 



1/2 



{asds + (TNdN) ■ (70) 



In the case of equal conductivities, as — a^, the effective 
thickness is simply the geometrical one: c?cff = ds + dN. 
This case corresponds to a uniform density of the super- 
conducting electrons, ns = tin, which implies a contin- 
uous distribution of the supercurrents, centered at the 
middle of the bilayer [this can be also seen from Eq. ( |6^ ) 
which yields zq — {ds+d^) /2 in the case as = aj^] . How- 
ever, in a more subtle situation when the conductivities 
are different, the density of the supercurrent experiences 
a jump at the SN interface; this nontrivial supercurrent 
distribution results in the nonequivalence of c?off to the 
geometrical thickness of the bilayer. 



VII. PERPENDICULAR UPPER CRITICAL 
FIELD 

Now we turn to calculating the upper critical field Hc2 
perpendicular to the plane of the bilayer. 

As in the case of the parallel critical field, we start with 
discussing the supercurrent distribution, which is now a 
function of the sample boundaries in the a;?/-plane, per- 
pendicular to the magnetic field H (the magnetic field 
is directed along the z-axis). The infinite bilayer under 
consideration can be thought of as a disk of a large ra- 
dius; let us assume a; = 0, ?/ = at the axis of the disk. 
Then the supercurrent distribution is axially symmetric, 
and, with the gauge chosen as A = [Hr] /2, the super- 
conducting phase must be constant, <^ = 0, which yields 
a simple result for the supercurrent velocity: v = eA/m. 

Near Hc2^ the superconducting correlations are small, 
\6\ <^ 1, and the Usadel equations can be linearized: 



y (-iV + 2eA)^ 



iE9 + A = 0, 
Ave = 0. 



(71a) 
(71b) 



The second of these equations is trivially satisfied because 
9{r) is axially symmetric. 

Th us, the Usadel equations reduce to the single 
Eq. (71a) for the pairing angle 9{r,E). Introducing 



11 



the cylindrical coordinates r ^ {z,p) and denoting 
P = —iVp + 2eA(p), we rewrite this equation as 



2 9z2 
Ds d^Os Ds 
2 dz^ 2 



-^P^On + iEeN ^0, (72a) 



+ iEOs + A = 0. 



(72b) 



We cannot solve these equations straightforwardly be- 
cause near the upper critical field, the order parameter 
A(p) is a nontrivial unknown function of the in-plane 
coordinate p (while the ^-dependence is absent due to 
the small thickness of the bilayer). In this situation, we 
employ the following approach. 

Averaging each of Eqs. ( [T^ over the thickness of the 
corresponding layer, we obtain 



d0N 




2dN 


dz 


z=Q 


Dn 


d9s 




2ds 


dz 


z=Q 


Ds 



IEOn 
- iEOs + A 



(73a) 
(73b) 



The averaged pairing angles entering the r.h.s. 
Eqs. (|7|) are 



of 



{P.E) 



rds 



dz 9n(z, p, E), 



(74a) 



esip, E) = ^ r dz es{z, p, E). (74b) 
ds Jo 



Substituting Eqs. (|7^) into the boundary conditions ( |2C 
(which should be linearized), we obtain a system of two 
differential equations for the function 9{p,E): 



TN 



D 



iE6 



N 



TS 



= OS 



iE9,, 



(75) 



From the vicinity of the superconductive transition it 
follows that the pairing angle 6 depends on the order 
parameter A linearly: 



0n{p,E) 
Os{p,E) 



A(p) 

ajv (E) ' 
A(p) 

as [EY 



(76a) 
(76b) 



where the fimctions aN(E) and as{E) are spatially in- 
dependent. Then Eqs. ( (Tq ) can be rewritten as 



^P^A(p) = 



Tjv V [Ej 



A(p), (77a) 



^P^A(p) = 



'^E+^f^-l 
Ts V "Jv [E) 



as{E) 



A(p). 
(77b) 



We see that the order parameter must be an eigenfunc- 
tion of the differential operator P^. Moreover, in order 
to obtain the largest value of _ffc2, we should choose the 
eigenfunction corresponding to the lowest eigenvalue (in 
complete analogy with Refs. [ 16 1^). The solution of 
the emerging eigenvalue problem is readily found thanks 
to its formal equivalence to the problem of determining 
the Landau levels of a two-dimensional particle with the 
"mass" 1/ D and the charge — 2e in the uniform magnetic 
field H directed along the third dimension. The lowest 
Landau level is DeH; the function as{E) is straightfor- 
wardly determined, 

(j^^ n u m{DNeH~iE) 

as{E) = DseH -lE + — -— — — -, 

TS [1 + TAT [Dn eH - iE)\ 

(78) 

and we substitute 9s (p, E) into the self-consistency equa- 
tion (^). The order parameter A(p) cancels out, and the 
resulting equation, which determines Hc2{T), can be cast 
into the form 



In- 



rpBCS 
T 



I \tstni^dJ \2 



Ts + tn 

1 



£s-S 



N 



1 - 



{£s-£n) +4/r5TAr 
Ss — 

yJ{£s-SNf + 4/TSTN, 



1 



1 



Ss + S 



N 



(Ss - SnY 



TSTN 



1 



2 AttT 



£s + £n — 



(Ss - En? + 

tstn 



(79) 
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where 





FIG. 9. Perpendicular upper critical field Hc2, normalized 
by its BCS value, versus pint at zero temperature. The up- 
per and lower graphs differ only in the scaling of the ordinate 
axis (normal and logarithmic, respectively). The inset shows 
Hc2ipint) on a wider scale over pint. 



£n = Dn eHc2 H 



(80) 



are iJ-dependent energies. The logarithmic term in the 
r.h.s. of Eq. ( [79| ) takes account of the finitencss of the 
Debye energy ujd, it becomes important only in the 
limit of a perfect interface (the Cooper limit), i.e., when 
tstni^d/ (rs + Tat) < 1. 

In the limit— Pint ^ oo, Eq.|-(£9|) yields the classical 
result of Makilla and de GennesliJfor the BCS case. 



In- 



T, 



BCS 



T 



2ttT 



(81) 



which is valid for bulk superconductors and superconduc- 
tive layers of arbitrary thickness (when the magnetic field 
is directed perpendicularly to them). Another immedi- 
ate consequence of Eq. (|79| ) is the critical temperature 
of the bilayer Tc, which can be found from the condition 
Hc2{Tc) = 0: in this case, Eq. (|7|) reproduces Eq. 

Equation ( [79| ) can be solved numerically at arbitrary 
values of the temperature T and the interface resistance 
Pint; the results for Hc2 are presented in Figs. ^, |l^. 

The analytical solution of Eq. ( [79| ) at zero temperature 
in the Anderson limit is presented below. 



A. Hc2 at zero temperature in the Anderson limit 




Pi„, = 80 and BCS curve 
Pm,= 125 



0.2 



0.4 



0.6 



T/T. 



FIG. 10. Temperature dependence of the perpendicular up- 
per critical field Hc2 at pint = 80, 110, 125, and in the BCS 
case. In each case, the critical field is normalized by its 
zero-temperature value, and the temperature is normalized 
by the corresponding T^. According to the results of Sec. Q 
the curves in the BCS and Anderson (pint = 80) limits coin- 
cide. At intermediate values of pint, the curves can lie both 
above (pint = HO) and below (pint = 125) the BCS curve. 



In the zero-temperature Anderson limit (defined by the 
conditions Ds eHc2 < l/rg, ZJ^v ei/c2 < I/tat), Eq. (|7£ 
yields 



(rg + Tjv) Dg 

H^cs - rsDs + tnD 



N 



^BCS L , f TsTNi^D 



2uj 



D 



\TS +TN 



(82) 



where the zero-temperature BCS value of the upper crit- 
ical field, as follows from Eq. (|8l|), is 



H. 



BCS 



A 



BCS 



$0 



2eDs 2t:^%cs 



(83) 



It is instructive to rewrite the perpendicular upper crit- 
ical field of the bilayer (B^) in the standard BCS-like form 



(84) 



where ^ is the bilayer's correlation length given by 
Eq. ( |69[ ) [the physical interpretation of this result for ^ 
precedes Eq. (|69|)]. 
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VIII. SNS, NSN, SNINS, NSISN, AND 
SUPERLATTICES 



ing discussion, it is convenient to write the condition of 
the small thickness of the bilayer as 



Our results for A, ns, Tc, Eg, v{E), and Hc2 (i.e., all 
the results except He) can be directly applied to more 
complicated structures such as SNS and NSN trilayers, 
SNINS and NSISN systems, and SN superlattices. 

Let us consider, for example, a symmetric S'A^S' trilayer 
consisting of two identical S layers of thickness ds sepa- 
rated by a iV layer of thickness 2(ijv. The SN interfaces 
can have arbitrary (but equal) resistances. As before, 
the z-axis is perpendicular to the plane of the structure. 
This trilayer can be imagined as composed of two iden- 
tical bilayers perfectly joined together along the sides. 
Indeed, the pairing angle 9 has zero z-derivative on the 
outer surfaces of the bilayers, thus producing the correct 
(symmetric in the z-direction) solution for 6 in the result- 
ing trilayer. Consequently, the symmetric SNS trilayer 
has exactly the same physical properties [A, rtg, Tc, Eg, 
v{E), Hc2\ as the S'iV bilayer considered in the present 
paper. The only point where the above reasoning fails 
is the calculation of the parallel critical field H^- In this 
case, the combination of the supercurrent distributions in 
the two bilayers does not yield the correct distribution in 
the resulting SNS trilayer, which implies that the Usadel 
equations for the two systems are different. 

Evidently, the above reasoning, based on the formal 
equivalence of the outer-surface boundary condition for 
the bilayer to the symmetry-caused condition in the mid- 
dle of the S'AfS' trilayer, also holds for symmetric NSN tri- 
layers (A^ layers of thickness c?Ar, 5 layer of thickness 2ds, 
identical SN interfaces) and SN superlattices (A^ layers 
of thickness 2dN, 5" layers of thickness 2ds, identical SN 
interfaces). Moreover, the same applies to systems com- 
posed of two bilayers in nonideal contact with each other: 
SNINS and NSISN (where / stands for an arbitrary po- 
tential barrier), because the presence of a potential bar- 
rier does not violate the applicability of the symmetry 
argument. Thus, all the results obtained for the bilayer 
(except He) are also valid for these structures. 



IX. DISCUSSION 

An essential property of the bilayer used throughout 
the paper is its small thickness. Now we shall argue that 
the bilayer studied in the experiment by Kasumov et alu 
(and to which our numerical results refer) can be con- 
sidered thin. The Usadel equations (|l^) imply that the 
characteristic space scale of the bilayer's properties vari- 
ation is y/DN,s / Eo for the N and S layers, respectively. 
However, the correct determination of the characteris- 
tic energy scale Eq is a nontrivial problem. Our results 
suggest that Eq is always smaller than the order param- 
eter A: in the BCS limit {pint oo), Eq approaches A, 
whereas in the opposite (Anderson) limit, Eq is deter- 
mined by the minigap Eg [see Eq. (34)]. For the foUow- 



dN.,s < 




BCS 



(85) 



The individual layers' thicknesses are d^ = 100 nm and 
ds = 5 nm. The third multiplier (the BCS correlation 
length) in the r.h.s. of the condition ( p5| ) equals 194 nm 
and 16 nm for the N and S layers, respectively. At the 
same time, each of the first two multipliers in the r.h.s. of 
the condition (85) exceeds unity. We can thus conclude 
that the bilayer can indeed be considered thin. 

Now we turn to a possible experimental application 
of our results. Our results provide a method for deter- 
mining Pint, a very important parameter of the bilayer 
which is nojt-directly measurable. By analyzing the ex- 
perimentalBEl values = 0.4 K and He ^ 0.1 T, we get 
Pint ~ 111 and Pint ~ 121, respectively. Within the ex- 
perimental accuracy of the bilayer's parameters, the two 
estimates for pint should be considered close. Interest- 
ingly, the value pint ~ 121 extracted from the measured 
value of He corresponds to the extremely narrow region 
of the steep drop in He{pint) (see Fig. ||). 

Finally, we wish to remark on a peculiarity of real 
systems which can be relevant when one compares our 
findings with an experiment. The point is that dur- 
ing the fabrication of a bilayer, the interface between 
S and A'^ materials cannot be made ideally uniform. In 
other words, the local interface resistance possesses spa- 
tial fluctuations. At the same time, as we have shown, 
the bilayer's properties are highly sensitive to the inter- 
face quality, which could lead to complicated behavior 
not reducing to the simple averaging of the interface re- 
sistance embodied in pint. One possibility could be a 
percolation- like proximity effect. We leave the study of 
inhomogeneity effects for further investigation. 



X. CONCLUSIONS 

We have studied, both analytically and numerically, 
the proximity effect in a thin S'A^ bilayer in the dirty 
limit. The layers were supposed to be thin enough to 
ensure uniform properties of each layer across its thick- 
ness. The strength of the proximity effect is governed by 
Pint , the resistance of the SN interface per channel. 

The quantities calculated were A, the order parameter; 
ns, the density of the superconducting electrons in the 
S layer; Te, the critical temperature; Eg and i^{E), the 
minigap in the density of states and the DOS itself; He 
and He2 , the critical magnetic field parallel to the bilayer 
and the upper critical field perpendicular to the bilayer. 

These quantities were calculated numerically over the 
entire range of pint. For this purpose, the characteristics 
of the bilayer were assumed tp be the same as in the ex- 
periment by Kasumov et aln that originally stimulated 
our research (Ta/Au bilayer, ds/^AT = 1/20). Inthehmit 
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of an opaque interface, A, ng, T^, iJc, and Hc2 approach 
their BCS values. At the same time, Eg does not coincide 
with the order parameter A, and Eg ^ when pi^t — > oo, 
although in general, the energy dependence of the DOS 
in the S and layers, vs{E) and vn{E), approaches the 
BCS and normal-metal results, respectively. 

The minigap Eg demonstrates nonmonotonic behav- 
ior as a function of pint- Analytical results for the two 
limiting cases of small and large pint show that in the An- 
derson limit, Eg increases with increasing pint, whereas in 
the limit of an opaque interface, Eg tends to zero. Thus, 
Eg reaches its maximum in the region of intermediate 

Pint- 

The most interesting case of relatively low interface 
resistance (the Anderson limit) has been considered an- 
alytically. The simple BCS relations between A, ns, Tc, 
He, Hc2 are substituted by similar ones with Eg stand- 
ing instead of A. The relation between the minigap Eg 
and the order parameter A in this limit is expressed by 
Eq. (34), implying that in the case where ts < tn, the 
BCS relations are strongly violated (by more than the 
order of magnitude for the above-mentioned Ta/Au bi- 
layer). The DOS in the S and A'^ layers coincide, show- 
ing BCS-like behavior with the standard peculiarity at 
E = Eg. It should be emphasized that absolute values of 
Pint corresponding to the Anderson limit can be large; for 
the Ta/Au bilayer this limit is already valid at pi^t < 80. 

All the results (except He) obtained for the bilayer 
also apply to more complicated structures such as SNS 
and NSN trilayers, SNINS and NSISN systems, and SN 
superlattices. 
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